Dynamics of metallic stripes in cuprates 
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We study the dynamics of metallic vertical stripes in cuprates within the three-band Hubbard 
model based on a recently developed time-dependent Gutzwiller approximation. As doping increases 
the optical conductivity shows transfer of spectral weight from the charge transfer band towards 
i) an incoherent band centered at 1.3eV, ii) a Drude peak, mainly due to motion along the stripe, 
Hi) a low-energy collective mode which softens with doping and merges with ii) at optimum doping 
in good agreement with experiment. The softening is related to the quasidegeneracy between Cu 
centered and O centered mean-field stripe solutions close to optimal doping. 
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The doping-dependent evolution from insulating be- 
havior to a strange metal in the superconducting cuprates 
emerges dramatically in the normal state optical conduc- 
tivity. Slightly doped cuprates show a small or no Drude 
peak and doping-induccd transfer of spectral weight from 
the charge-transfer (^2eV) to a mid-IR (MIR) band at 
'^O.SeVllI, ||. Upon further doping the system progres- 
sively metallizes as is evident from the prominent Drude- 
like peak that develops at zero energy. A remarkable 
effect of doping is that the MIR band strongly softens 
and merges with the Drude peak, resulting in a feature 
that cannot be fitted by a conventional Drude model. 
A variety of alternative theories ||, ^ have been pro- 
posed in order to describe this feature. Clearly the iden- 
tification of this low-energy MIR (LEMIR) band is of 
paramount importance to understand the physics of these 
materials. The low-doping behavior has been explained 
in terms of the random-phase-approximation (RPA) elec- 
tronic excitations of single-hole Hartree-Fock states in 
Cu02 layers[^, but the moderate doping behavior (the 
softening of the LEMIR band) could not be explained 
due to difficulties with the HF ground state. 

The softening of the LEMIR band in La2-a;Sra;Cu04 
(LSCO) is accompanied by the appearance of another 
(much less discussed) band at 1.3eV||l|, H, [g). This high 
energy MIR (HEMIR) band is well pronounced in op- 
tical absorption through LSCO thin films [||, and elec- 
tron energy loss spectroscopy where it develops as a 
function of doping. Moreover LEMIR and HEMIR are 
also detected in photodoped experiments on LSCO0. 
The HEMIR has not been clearly resolved by reflectivity 
in YBa2Cu306+5 (YBCO), but a strong broad feature 
at the right energy appears in photodoped transmission 
experiments As far as we know no microscopic expla- 
nation of the HEMIR exists so far. 

Another important aspect of layered cuprates that has 
emerged in the last years is the rearrangement of doped 
holes in antiferromagnetic (AF) domain walls ^ 
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FIG. 1: Charge density minus undoped charge density as a 
function of atomic position in the direction perpendicular to 
the d = 4 stripe for the BC solution at a; = 0.17 (a) and the 
SC solution at a; = 0.22 (c). (b) shows the transition charge 
for the LEMIR collective mode of the panel (a) solution. 



O, 01 H. These one-dimensional (ID) structures called 



stripes were predicted by mean- field theories 1 15 



Recently we have presented a computation of metallic 
mean- field stripes]!^ within an unrestricted Gutzwiller 
approximation ||l7[ (GA). The behavior of the magnetic 
incommensurability e = l/(2d)||ll[ |l2[ ^ (d is the 
distance between charged stripes in units of the lat- 
tice constant), chemical potential [l^, and transport 
experiments ||20|, ^ as a function of doping have been 
explained in a parameter free wayp6|. 

In this work we investigate the optical conductivity 
in terms of linear excitations around the metallic mean- 
field stripes of Ref. [Q. RPA fluctuations are added 
through a recently developed time-dependent GA la- 
beled GA-l-RPAp^. The stripes show a strong Drude 
component in the stripe direction, of magnitude similar 
to the experimental one. In addition the single-holc[^ 
MIR band splits into two bands in good accord with the 
HEMIR and LEMIR (Fig. |). The HEMIR band can be 
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understood in terms of incoherent transitions between 
band states of a stripe at mean-field level (Fig. ^). The 
LEMIR band is a collective mode associated with soft lat- 
eral displacements of the stripe. As doping increases the 
collective mode shifts to lower energy and merges with 
the Drude part, providing an explanation for the soften- 
ing of the LEMIR band in cuprates. This softening is 
due to the quasidegeneracy between 0-centered and Cu- 
centered stripes found in many approaches Q and 

the small energy barrier between these solutions (Fig. |^). 
The soft collective mode is a good candidate to mediate 
pairing between holes and to produce non-Fermi liquid 
anomalies in the normal state. A related scenario was 
presented in Ref. ^ but the origin of the softening was 
not clear. 

We use a three-band model for cuprates with the same 
LDA parameter set[^ as in Ref. [ p^ . 

We start by reviewing some mean- field results |, 
p6t in order to illustrate the doping dependent charge 
and spin textures on which the following computations 
are based. At very dilute doping {x ^ 0.03) due to the 
long-range Coulomb interaction (not included in our cal- 
culations) each hole will be close to an acceptor prevent- 
ing the formation of stripes. The lowest energy one-hole 
solution consists of a self-trapped state similar to the 
Zhang- Rice state as found in HF|^, |2^ . 

As doping increases, the donor potential becomes more 
uniform and screened, favoring the formation of stripes. 
Experiment shows that stripes are parallel to the Cu- 
O bond except at dopings 0.03 < x 
agonal stripes have been observed 
an intermediate state between the isolated polarons and 
the vertical stripes and probably also require the long- 
range Coulomb interaction to be stable. For simplicity 
we skip this phase and consider vertical metallic solu- 
tions. Weak one-dimensional instabilities probably rele- 
vant at low temperatures are intentionally suppressed in 
our clusters due to finite size effects. 

In Ref. we have shown that the most favorable low- 
doping metallic mean-field stripe is centered on O sites 
bridging two vertical Cu legs, denoted as bond-centered 
(BC) in one band models|^, Fig. ^ (a) shows a 

cross section of the charge modulation perpendicular to 
the stripe. For larger concentrations the BC solution 
becomes degenerate with the site centered (SC) one at 
xo ~ 0.21 whose charge profile is shown in Fig. |^(c). 

The mean-field bands for a BC stripe lattice as re- 
ported in Ref. [|l6| are shown in Fig. ||. Here we discuss 
in more detail the symmetry which determines the selec- 
tion rules for optical transitions. Roughly speaking the 
flat bands labeled S and P correspond to symmetric (S) 
and antisymmetric (P) combinations of orbitals centered 
on the two legs of Cu that form the core of the stripe 
(sites 2 and 3 in Fig. |l|). The band crossing the chemi- 
cal potential (hereafter the "active band") is due to the 
orbitals centered on the core O leg of the stripe [at 2.5 
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FIG. 2; Electron mean-field bands measured from the chemi- 
cal potential for the momentum in the direction of the stripe 
{d = 7, X = 0.071). The arrows are the lowest energy dipole 
allowed mean-field transitions in a d x 4 cluster labeled by 
energy and polarization. Notice that the x polarization is 
perpendicular to the stripe. We also plot the energies of the 
AF insulating bands at momentum k — (0, vr) (full dots) mea- 
sured from the same reference energy 



in Fig. |^(a)]. The antisymmetric Cu orbital combination 
mixes with the core O orbital pushing upwards (down- 
wards) the P band (active band) close to the edge of 
the stoichiometric bands (marked by the full dots). The 
lower bands are of mainly O character. All optical tran- 
sitions in the x direction indicated in Fig. ^ are of course 
between even and odd states with respect to the stripe 
central axis. 

Dynamical properties are computed within a real space 
implementation of the GA-I-RPA method [p2[. In this 
context it is essential to use GA-t-RPA rather than con- 
ventional HF-I-RPA theory due to the above-mentioned 
fact that the starting point mean-field solution matches 
the experimentally observed doping-induccd incommen- 
surability. As a bonus GA+RPA has been shown to be 
more accurate than HF+RPAp^. Note that stripes in 
LSCO show a dynamic character on a scale of a few 
meV[|l^ which is not capture by our starting-point mean- 
field ground state. We expect this to affect the spectra 
at very low energies as discussed in Ref. but not on 
the scale of the transitions reported Fig. ||. 

Fig. ^ reports the optical conductivity within 
GA-I-RPA for various dopings. The RPA optical con- 
ductivity for the single-hole solution appropriate at low 
doping, as discussed above, reproduces the results of 
HF-|-RPA|^: Formation of a doping induced MIR band 
close to 0.5eV and doping induced transfer of spectral 
weight from the charge transfer band to the MIR region 
in agreement with experiment in this doping range Q]. 

For distant stripes {d — 7) the single-hole MIR band 
now splits into two bands. The one at higher energy 
is a band of incoherent particle-hole excitations close to 
1.3eV which provides a theoretical explanation for the 
HEMIR. The position of this band is nearly indepen- 
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FIG. 3: Optical conductivity labeled by doping, system size 
and, in the case of stripes, interstripe distance. The units of 
conductivity are given by ctq = 3.6 x 10^ (ficm)"^ with a back- 
ground dielectric constant £b = 2.6 (see Ref. Q and caption 
of Fig. P). The curve labeled x = 0.028 corresponds to the 
single-hole solution. For larger dopings the figure is an aver- 
age over the electric field directions parallel and perpendicular 
to the BC stripes. The inset shows the low energy spectra ex- 
cluding the Drude component. We used a Lorentzian broad- 
ening of 0.2eV. 



dent of doping and can be understood in terms of transi- 
tions within the stripe band structure at mean-field level 
shown in Fig. ^. Indeed the HEMIR is mainly formed by 
the 1.33eV(x) and 1.34eV(?/) mean field transitions with 
similar oscillation strengths. This negligible renormal- 
ization of the mean field transitions by the RPA is char- 
acteristic for incoherent particle-hole excitations. The 
other MIR band, shown also in the inset of Fig. ||, is a 
low-energy collective mode and has no mean-field coun- 
terpart. In fact the 0.85eV(a;) transition which has a 
strong oscillator strength in mean-field does not show up 
in RPA and instead the low-energy collective mode ap- 
pears which is also polarized perpendicular to the stripes. 
A similar mode was found in a study of stripes within the 
t — t' — t" — J model . The other transitions reported 
in Fig. ^ have much smaller spectral weights. 

To characterize the LEMIR band we compute the 
so called transition charges nf = (0|ni|//) where la- 
bels RPA excitations, nj' is proportional to the time- 
dependent charge fluctuation Sni that would occur at 
frequency lo^ if the state v were weakly excited 

In Fig. |^(b) we show the charge fluctuation associated 
with the LEMIR mode. It is very similar to the differ- 
ence in charge density between the BC and SC stripe, 
indicating that it corresponds to lateral displacements of 
the stripe. In fact, if we approximate the charge modula- 
tion by cos(q'cdw'' + d) with gcdw = 27r(2e, 0), this oscilla- 
tion can be interpreted as a time-dependent fluctuation 
of the phase and thus the LEMIR excitation can be 
identified with a phason. Optically active phasons have 
zero momentum, but naturally a band of phasons ex- 
ist with a well-defined dispersion relation. In continuum 



FIG. 4: Mean field energy per stripe cell (8 Cu's) as a function 
of the collective phase 9 of the CDW for various doping levels 
and for d = 4 stripes in a 4d x 12 system. 

models phasons are massless Goldstone modes whereas 
here the commensurability of Qcdw with the lattice makes 
them have a finite energy at zero momentum. This en- 
ergy however is small and decreases as doping increases. 
Since the BC and SC state become quasidegenerate at 
xo it is natural to expect that the phason softening is 
related to the quasidegeneracy between these states. 

In order to substantiate this idea we compute the en- 
ergy landscape for "intermediate" solutions constrained 
to be stripes centered between Cu and O. For each in- 
termediate state we perform a Fourier analysis and ex- 
tract the phase 9 of the first Fourier component of the 
charge-density-wave (CDW) modulation. Fig. ^ shows 
the energy for d = 4 stripes and different dopings. This 
provides an upper bound for the energy along the path 
connecting BC {9 = 7r/4) and SC {9 = 0) solutions where 
9 plays the role of a collective coordinate. The curves are 
periodic in 9 with period 27r/d corresponding to a transla- 
tion by one elementary unit cell. Remarkably the curve 
acquires an extra periodicity close to optimum doping 
corresponding to the previously found quasidegeneracy 
between SC and BC solutions[^ || ||. 

RPA is essentially an harmonic approximation of the 
energy landscape around the mean field solution (SC or 
BC). The energy squared of each RPA mode is propor- 
tional to the curvature (or "stiffness" ) of the correspond- 
ing parabolic energy approximation when the system is 
displaced from the stationary state in the direction of the 
mode eigenvector (in our case parametrized by the collec- 
tive coordinate 9). Fig. ^ shows that as doping increases 
the stiffness decreases showing explicitly that the soften- 
ing of the LEMIR feature is due to the quasidegeneracy 
between BC and SC stripes. Of course the problem is 
very anharmonic close to and RPA provides only a 
rough estimate to the phason energy. Moreover, since 
the barrier is strongly reduced close to xq we expect that 
anharmonic corrections will make the phason even softer 
and the "true" ground state will be a fluctuating mixture 
of BC and SC solutions. 
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FIG. 5: Experimental[Q] and theoretical optical conductivity 
spectral weight integrated up to 1.5 eV vs. x. The spectral 
weight is converted to an effective number of electrons A'cft 
as in Ref. The background dielectric constant is our 
only free parameter and was adjusted to make the theoretical 
and experimental intensities to have an overall agreement. 
We also show the computed Drude and regular contributions 
in each direction. The dashed lines indicate the region of 
metastability of the d = 4 solution. The points not joined by 
lines correspond to the d = 3 solution which becomes more 
stable in that region. 

In Fig. ^ we show the optical conductivity spectral 
weight integrated up to an energy of 1.5 eV [TVos (1.5eV)] 
in comparison with experimental data from Ref. |l[] . The 
lower curves report the regular (w > 0) and Drude 
(w = 0) contributions to iVoff (1.5eV) in each direction. 

For X < 1/8 the number of stripes, which 
act as low-energy light absorbers, increases linearly 
with doping keeping its electronic structure practi- 
cally unchanged p6|. Accordingly A^eff(l-5eV) increases 
roughly linearly with doping. For 1/8 > x > x\ 
{x\ ~ 0.22 — 0.23) the number of absorbers (stripes) get 
locked and the evolution of the spectral weight is related 
to changes in the electronic structure of each stripe. Es- 
pecially the shift of the chemical potential from the cen- 
ter of the active band||l^ leads to a depression of Drude 
weight in the stripe direction (joined circles) and results 
in a slower increase of spectral weight with doping which 
correlates well with experiment. 

For X > x\ d — "i stripes become the mean-field ground 
state. These solutions tend to have a larger spectral 
weight due to the strong Drude weight in both directions 
(unjoined points in Fig. The nature of the ground 
state however is not clear in this overdoped regime and 
therefore our results become more qualitative than quan- 
titative. Indeed d = 3 stripes have not been observed in 
LSCO and theoretically a combination of different solu- 
tions would likely be more appropriate which is be- 
yond the scope of this work. 

To conclude we have computed the optical conduc- 
tivity of metallic stripes within a GA+RPA approach. 
Striped domain walls induce two MIR excitations: At 
around 1.3eV a HEMIR band appears which is related 



to interband transitions within the stripe band structure. 
Further on we found a collective mode (i.e. the LEMIR 
band) which softens as a function of doping due to the 
suppression of the energy barrier between quasidegener- 
ate BC and SC stripe solutions. These features are in 
good agreement with experiments for positions and rel- 
ative intensities with parameters fixed by first principle 
computations as in Ref. |l6j . It is worth speculating that 
scattering of holes moving along and perpendicular to the 
stripes with the soft collective mode may be responsible 
for both the anomalous normal state behavior and the 
superconducting pairing in cuprates. 

After this work was completed Ref. |^ was posted in 
which by studding temperature effects, MIR features in 
isostructural La2Ni04.i33 were also linked to stripes. 
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Forschungsgemeinschaft . 
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